A GENERALIZATION OF THE LOWNER-JOHN'S 
ELLIPSOID THEOREM 



JEAN B. LASSERRE 



Abstract. We address the following generalization P of the Lowner- 
John ellipsoid problem. Given a (non necessarily convex) compact set 
K C R" and an even integer d £ N, find an homogeneous polynomial g 
of degree d such that K C G := {x : g(x) < 1} and G has minimum 
volume among all such sets. We show that P is a convex optimization 
problem even if neither K nor G are convex! We next show that P has 
a unique optimal solution and a characterization with at most ("^^~^) 
contacts points in KnG is also provided. This is the analogue for d > 2 
of the Lowner- John's theorem in the quadratic case d — 2, but impor- 
tantly, we neither require the set K nor the sublevel set G to be convex. 
More generally, there is also an homogeneous polynomial g of even de- 
gree d and a point a £ K" such that K C Ga {x : g{x — a) < 1} 
and Ga has minimum volume among all such sets (but uniqueness is 
not guaranteed). Finally, we also outline a numerical scheme to approx- 
imate as closely as desired the optimal value and an optimal solution. 
It consists of solving a hierarchy of convex optimization problems with 
strictly convex objective function and Linear Matrix Inequality (LMI) 
constraints. 



1. Introduction 

"Approximating" data by relatively simple geometrical objects is a fun- 
damental problem with many important applications and the ellipsoid of 
minimum volume is the most well-known of the associated computational 
techniques. 

Indeed, in addition to its nice properties from the viewpoint of appli- 
cations, the ellipsoid of minimum volume is also very interesting from a 
mathematical viewpoint. Indeed, if K C M" is a convex body, computing 
an ellipsoid of minimum volume that contains K is a classical and infamous 
problem which has a unique optimal solution called the Lowner- John's el- 
lipsoid. In addition, John's theorem states that the optimal ellipsoid Vt 
is characterized by s contacts points € K n fi, and positive scalars Aj, 
i = 1, . . . , s, where s is bounded above by n{n -|- 3)/2 in the general case and 
s < n{n + l)/2 when K is symmetric; see e.g. Ball OHj. In particular, and 
in contrast to other approximation techniques, computing the ellipsoid of 
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minimum volume is a convex optimization problem for which efficient tech- 
niques are available; see e.g. Calafiore [7] and Sun and Freund |32J for more 
details. For a nice recent historical survey on the Lowner-John's ellipsoid, 
the interested reader is referred to the recent paper by Henk [l3j and the 
many references therein. 

As underlined in Calafiore [7], "The problem of approximating observed 
data with simple geometric primitives is, since the time of Gauss, of fun- 
damental importance in many scientific endeavors". For practical purposes 
and numerical efficiency the most commonly used are polyhedra and ellip- 
soids and such techniques are ubiquitous in several different area, control, 
statistics, computer graphics, computer vision, to mention a few. For in- 
stance: 

• In robust linear control, one is interested in outer or inner approxima- 
tions of the stability region associated with a linear dynamical system, that 
is, the set of initial states from which the system can be stabilized by some 
control policy. Typically, the stability region which can be formulated as 
a semi-algebraic set in the space of coefficients of the characteristic poly- 
nomial, is non convex. By using the Hermite stability criterion, it can be 
described by a parametrized polynomial matrix inequality where the pa- 
rameters account for uncertainties and the variables are the controller co- 
efficients. Convex inner approximations of the stability region have been 
proposed in form of polytopes in Nurges [23], ellipsoids in Henrion et al. 
|15) . and more general convex sets defined by Linear Matrix Inequalities 
(LMIs) in Henrion et al. [l7], and Karimi et al. |18] . 

• In statistics one is interested in the ellipsoid ^ of minimum volume cover- 
ing some given k oi m data points because ^ has some interesting statistical 
properties such as affine equivariance and positive breakdown properties [9] . 
In this context the center of the ellipsoid is called the minimum volume el- 
lipsoid (MVE) location estimator and the associated matrix associated with 
^ is called the MVE scatter estimator; see e.g. Rousseeuw [31] and Croux et 
al. [9]. 

• In pattern separation, minimum volume ellipsoids are used for separating 
two sets of data points. For computing such ellipsoids, convex programming 
techniques have been used in the early work of Rosen [28] and more modern 
semidefinite programming techniques in Vandenberghe and Boyd [34] . Sim- 
ilarly, in robust statistics and data mining the ellipsoid of minimum volume 
covering a finite set of data points identifies outliers as the points on its 
boundary; see e.g. Rousseeuw and Leroy [31j. Moreover, this ellipsoid tech- 
nique is also scale invariant, a highly desirable property in data mining which 
is not enjoyed by other clustering methods based on various distances; see 
the discussion in Calafiore [3, Sun and Freund [32] and references therein. 

• Other clustering techniques in computer graphics, computer vision and 
pattern recognition, use various (geometric or algebraic) distances (e.g. the 
equation error) and compute the best ellipsoid by minimizing an associated 
non linear least squares criterion (whence the name "least squares fitting 
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ellipsoid" methods). For instance, such techniques have been proposed in 
computer graphics and computer vision by Bookstein [6] and Pratt [25], in 
pattern recognition by Rosin [29], Rosin and West [30], Taubin [33], and in 
another context by Chernousko [8]. When using an algebraic distance (like 
e.g. the equation error) the geometric interpretation is not clear and the re- 
sulting ellipsoid may not be satisfactory; see e.g. an illuminating discussion 
in Gander et al. [12]. Moreover, in general the resulting optimization prob- 
lem is not convex and convergence to a global minimizer is not guaranteed. 

So optimal data fitting using an ellipsoid of minimum volume is not only 
satisfactory from the viewpoint of applications but is also satisfactory from 
a mathematical viewpoint as it reduces to a (often tractable) convex opti- 
mization problem with a unique solution having a nice characterization in 
term of contact points in K fl fi. In fact, reduction to solving a convex 
optimization problem with a unique optimal solution, is a highly desirable 
property of any data fitting technique! 

A more general optimal data fitting problem. In the Lowner-John 
problem one restricts to convex bodies K because for a non convex set K 
the optimal ellipsoid is also solution to the problem where K is replaced 
with its convex hull co(K). However, if one considers sets that are more 
general than ellipsoids, an optimal solution for K is not necessarily the 
same as for co(K), and indeed, in some applications one is interested in 
approximating as closely as desired a non convex set K. In this case a non 
convex approximation is sometimes highly desirable as more efficient. 

For instance, in the robust control problem already alluded to above, in 
Henrion and Lasserre [16] we have provided an inner approximation of the 
stability region K by the sublevel set G = {x : ^(x) < 0} of a non convex 
polynomial g. By allowing the degree of g to increase one obtains the con- 
vergence vol(G) vol(K) which is impossible with the convex polytopes, 
ellipsoids and LMI approximations described in [23 1 ITS j [T7 1 [T8]. 

So if one considers the more general data fitting problem where K and /or 
the (outer) approximating set are allowed to be non convex, can we still infer 
interesting conclusions as for the Lowner-John problem? Can we also derive 
a practical algorithm for computing an optimal solution? If the answer is 
negative, such a generalization would then be useless or sterile (or at least 
questionable) . 

The purpose of this paper is to provide results in this direction that can 
be seen as a non convex generalization of the Lowner- John's problem but, 
surprisingly, still reduces to solving a convex optimization problem with a 
unique solution. 

Some works have considered generalizations of the Lowner-John's prob- 
lem. Fo instance, Giannopoulos et al. [lOj have extended John's theorem for 
couples (Ki,K2) of convex bodies when Ki is in maximal volume position 



4 



JEAN B. LASSERRE 



of Ki inside K2, whereas Bastero and Romance [5] refined this result by 
allowing Ki to be non-convex. 

In this paper we consider a different non convex generalization of the 
Lowner-John ellipsoid problem, with a more algebraic flavor. Namely, we 
address the following two problems Pq and P. 

Pq.' Let K C M" be a compact set (not necessarily convex) and let d he 
an even integer. Find an homogeneous polynomial g of degree d such that 
its suhlevel set Gi := {x : ^(x) < 1} contains K and has minimum volume 
among all such sub level sets with this inclusion property. 

P; Let K C M" be a compact set (not necessarily convex) and let d he an 
even integer. Find an homogeneous polynomial g of degree d and a € M" 
such that the suhlevel set G\ := {x : (/(x — a) < 1} contains K and has 
minimum volume among all such sub level sets with this inclusion property. 

Necessarily (7 is a nonnegative homogeneous polynomial otherwise Gi 
and Gf are not bounded. Of course, when d = 2 then g is convex (i.e., Gi 
and Gf are ellipsoids) because every nonnegative quadratic form defines a 
convex function, and g is an optimal solution for problem P with K or its 
convex hull co(K). That is, one retrieves the Lowner- John's problem. But 
when d > 2 then Gi and Gf are not necessarily convex. For instance, take 
K = {x : ^(x) < 1} where g is some nonnegative homogeneous polynomial 
such that K is compact but non convex. Then g is an optimal solution for 
problem Pq with K and cannot be optimal for co(K); a two-dimensional 
example is (x, y) i-)- g{x, y) := x^y"^ + e{x^ + y^) for e > sufficiently small. 

Contribution. Our contribution is to show that problem Pq and P are 
indeed natural generalizations of the Lowner-John ellipsoid problem in the 
sense that: 

- Pq also has a unique solution g* . 

- A characterization of g* also involves s contact points in K n Gi, where 
s is now bounded by ("'^^~^)- 

And so when d = 2 we retrieve the symmetric Lowner-John problem as a 
particular case. In fact it is shown that Pq is a convex optimization problem 
no matter if neither K nor Gi are convex! 

We use an intermediate and crucial result of independent interest. Namely, 
the Lebesgue-volume function g 1— > f{g) := vol(Gi) is a strictly convex 
function of the coefficients of g, which is far from being obvious from its 
definition. 

Concerning the more general problem P, we also show that there is an op- 
timal solution (g* ,a*) with again a characterization which involves s contact 
points in K n G^ , but uniqueness is not guaranteed. 

Again and importantly, in both problems Pq and P, neither K nor G^ 
are required to be convex. 
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On the computational side. Even though Pq is a convex optimization 
problem, it is hard to solve because even if K would be a finite set of of 
points (as is the case in statistics applications of the Lowner-John prob- 
lem) and in contrast to the quadratic case, evaluating the (strictly convex) 
objective function, its gradient and Hessian can be a challenging problem, 
especially if the dimension is larger than n = 3. So this is one price to 
pay for the generalization of the Lowner-John's ellipsoid problem. (Notice 
however that if K is not a finite set of points then even the Lowner-John's 
ellipsoid problem is also hard to solve because for more general sets K the 
inclusion constraint K C ^ (or conv(K) C ^) can be difficult to handle.) 

However, we can still approximate as closely as desired the objective func- 
tion as well as its gradient and Hessian by using the methodology developed 
in Henrion et al [U], especially when the dimension is small n = 2, 3 (which 
is the case in several applications in statistics). 

Moreover, if K is a (compact) basic semi- algebraic with an explicit de- 
scription {x : 5j(x) > 0, J = 1, . . . ,m} for some polynomials (gj) C IR[x], 
then we can use powerful positivity certificates from real algebraic geom- 
etry to handle the inclusion constraint Gi D K in the associated convex 
optimization problem. 

Therefore, in this context, we also provide a numerical scheme to approx- 
imate the optimal value and the unique optimal solution of Pq as closely as 
desired. It consists of solving a hierarchy of convex optimization problems 
where each problem in the hierarchy has a strictly convex objective function 
and a feasible set defined by Linear Matrix Inequalities (LMIs). 

2. Notation, definitions and preliminary results 

2.1. Notation and definitions. Let M[x] be the ring of polynomials in the 
variables x = (xi, . . . ,x„,) and let M[x]^ be the vector space of polynomials 
of degree at most d (whose dimension is s{d) := ("^'^))- For every d G N, 
let := {a G N" : \a\ (= J^iOi) < d}, and let v/x) = (x°), a G N", be 
the vector of monomials of the canonical basis (x") of M[x](i. 
A polynomial / G ^[x]^; is written 

x^/(x) = ^ /„x-, 

for some vector of coefficients f = (/„) G M^'^'^^. A polynomial / G M[x]rf is 
homogeneous of degree d if /(Ax) = A'^/(x) for all x G M" and all A G M. 

Let us denote by P[x]rf C M[x]rf, d G N, the set of nonnegative and 
homogeneous polynomials of degree d such that their sublevel set Gi := 
{x : ^(x) < 1} has finite Lebesgue volume (denoted vol(Gi)). Notice that 
necessarily d is even and P[x]rf. 

For d G N and a closed set K C M", denote by Crf(K) the convex cone 
of all polynomials of degree d that are nonnegative on K, and denote by 
A^(K) the Banach space of signed Borel measures with support contained 
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in K (equipped with the total variation norm). Let M(K) c M.(K) be the 
convex cone of finite Borel measures on K. 

In the EucUdean space R" we denote by (•, •) the usual duality bracket. 



2.2. Some preliminary results. We first have the following result: 



Proposition 2.1. The set P[yi\d is a convex cone. 



Proof. Let g,h with associated sublevel sets Gi = {x : ^(x) < 

1} and Hi = {x : /i(x) < 1}. For A G (0,1), consider the nonnegative 
homogeneous polynomial 9 := Xg+{l — \)h € M[x]rf, with associated sublevel 
set 



Gi := {x : 0(x) < 1} = {x : Xg{:x.) + (1 - A)/i(x) < 1}. 



Write Bi = U Gf where 9^ = Gi n {x : ^(x) > /i(x)} and G? = 
Gi n {x : (/(x) < /i(x)}. Observe that x € G} implies /i(x) < 1 and so 
G} C Hi. Similarly x G Gf implies g{x.) < 1 and so Gi C Gi. Therefore 
vol(Gi) < vol(Gi) + vol(Hi) < oo. And so 6* G P[x]d. □ 



With y G M and G M[x] let Gj^ := {x : c/(x) < y}. The following 
intermediate result is crucial and of independent interest. 



Lemma 2.2. Let g G P[x]d. Then for every y > 0: 



n/d r 

Proof. Asg e P[x]d, and using homogeneity, vol(Gi) < oo implies vol(Gj;) < 

oo for every y > 0. Let / : M — ?• M be the function y i— >■ f{y) := vol(Gy). 
Since g is nonnegative, the function / vanishes on (— oo,0]. Its Laplace 
transform : C — )■ C is the function 



/•oo 

X^Cf{X):= exp{-Xy)f{y)dy, 3?A > 0. 
Jo 
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Observe that whenever A G M with A > 0, 

£/(A) = / exp(-Ay) / dx] dy 

Jo \J{x:g{x)<y} / 

= / / ex.p{—Xy)dy dx [by Fubini's Theorem] 

\7g(x) J 

1 
A 



exp(— Ag(x)) dx 
- / exp{—g{X^^^^^x))dx [by homogeneity] 
J^^^exp(-5(z))dz [byA-VMx^z] 



A 

A1+ 



exp(-g(z)) dz 
T{l+n/d) ^W'*^^)- 
And so, by analyticity of the Laplace transform Cf, 

yU/d r- 

= rn _L u\ / exp(-5(x)) dx, y > 0, 
r(l + n/d) Jj^n 

which is the desired result. □ 

And we also conclude: 
Corollary 2.3. g G P[x]d <:=^ / exp{—g{x.))dx < oo. 

Proof. The implication ^ follows from Lemma 12.21 For the reverse impli- 
cation consider the function 

A h(X) := — / exp(-5(x)) dx, < A G M. 

Proceeding as in the proof of Lemma 12.21 but backward, it follows that h 
is the Laplace transform of the function y ^ vol(Gy) and the proof is 
complete. □ 

Formula ()2.ip relating the Lebesgue volume Gi with J" exp(— g') is already 
proved (with a different argument) in Morozov and Shakirov |201 [21] where 
the authors want to express the non Gaussian integral J exp(— in terms 
of algebraic invariants of g. 

Sensitivity analysis and convexity. We now investigate some properties 
the function / : P[x]rf — )• M defined by 

(2.2) g ^ fig) := vol(Gi) = [ dx, g e P(x]d, 

J{x:a(x)<l} 
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i.e., we now view vol(Gi) as a function of the vector g = (ga) € M^'-'^^ of 
coefficients of g in tlie canonical basis of homogeneous polynomials of degree 
d (and^(d) = C^^-^)). 

Theorem 2.4. The Lebesgue-volume function f : P(x]rf — ?> M defined in 
112. ^) is strictly convex with gradient V/ and Hessian V^/ given by: 



(2.3) 



dfig) 



r(i + n/d) L 



x" exp(-c/(x)) dx, 



for all a G N^, \a\ = d. 
(2.4) l^*"' - 



1 



r(l +n/d) 

for alia, p e NJJ, |a| = |/3| = d. 
Moreover, we also have 



exp(-g(x)) (ix, 



(2.5) 



5(x) exp(-5(x))dx 



n 



exp(-5f(x)) dx. 



Proof. By Lemma [2^2] /(^r) = r(l + n{d)) ^ f^„ exp(— 5f(x))(ix. Let p,q € 
P[x]rf and a E [0, 1]. By convexity of u i-7> exp(— u), 

/(a^f + (1 - a)g) < / [a exp(-fif(x)) + (1 - a) exp(-g(x)) ] dx 
= a/(5) + (l-a)/(g), 

and so / is convex. Next, in view of the strict convexity of n i-^ exp(— u), 
equality may occur only if ^(x) = g(x) almost everywhere, which implies 
g = q and which in turn implies strict convexity of /. 

To obtain (|2.3p - (|2.4p one takes partial derivatives under the integral sign, 
which in this context is allowed. Indeed, write g in the canonical b clSlS clS 
ai^) = Yl\a\=d9a^"- For every a G with \a\ = d, let e« = (e«(^)) G M^^'') 
be such that ea(/3) = 5/3=a (with 6 being the Kronecker symbol). Then for 
every t > 0, 



f{g + te^) - fig) 



exp{-g) 



exp(— tx") — 1 



V 



dx. 



Notice that for every x, by convexity of the function t i— )• exp(— tx"), 



liimh(t,x.) = inf^(t,x) = exp(— tx°)|,_n 
40 t>o l'-*^ 



-X 



because for every x, the function t i->- Tp{t,x) is nondecr easing; see e.g. 
Rockafellar [27^ Theorem 23.1]. Hence, the one-sided directional derivative 
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f'{g;ea) in the direction Ca satisfies 

f'{Q\en) = lim ~^ ^^"^ /(sO _ f gxn(— q) ■(/'(t, x) dx 

' m t 40 Ju" ^ ) 

= / exp(— 5^) lim^/;(t,x) (i^(x) = / — x" exp(— (7) (ix, 

JR" 40 JlRn 

where the third equahty fohows from the Extended Monotone Convergence 
Theorem |2l 1.6.7]. Indeed for all t < Iq with to sufficiently small, the 
function ^(t^ ■) is bounded above by V'l^O) ■) and /jg„ exp{—g)ip{to,x.)dfj, < 00. 
Similarly, for every t > 



S(t,x) 

and by convexity of the function t 1— > exp(tx") 



limf(t,x) = inff(t,x) = expftx'^ n.n 
40 t>o 1'-*^ 

Therefore, with exactly same arguments as before, 

.// X ,. fjg - teg) - /(g) 
fig; -eg) = lim 



and so 

df{g) 



/ x°'exp{-g)dx = -f'{g;ea), 
= - x" exp(-5r)dx. 



for every a with \a\ < d, which yields (j2.3p . Similar arguments can used for 
the Hessian V"^ f{g) which yields ()2.4p . 

To obtain (j2.5p observe that g 1— )• := J exp( — (7)(ix, g G P[x]rf, is a 
positively homogeneous function of degree — n/d, continuously differentiable. 
And so combining (j2.3p with Euler's identity yields 

/ exp(-g(x))(ix = --h{g) 
a j^n a 

= {'Vh{g),g) [by Euler's identity] 
g(x) exp(-g(x))(ix. 

□ 



Notice that convexity of / is not obvious at all from its definition (j2.2p 
whereas it becomes almost transparent when using formula (12. ip . 
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2.3. The dual cone of Cd(K). For a convex cone C C M", the convex cone 

C* := {y : (y,x) > Vx e C}, 

is called the dual cone of C, and if C is closed then (C*)* = C. 

Recall that for a set K C M", Cd(K) denotes the convex cone of polyno- 
mials of degree at most d which are nonnegative on K. We say that a vector 
y G M^(°'^ has a representing measure (or is a d-truncated moment sequence) 
if there exists a finite Borel measure (j) such that 

x" d(P, Vq e N^. 

We will need the following characterization the dual cone Cd(K.)* ■ 

Lemma 2.5. Let K C be compact. For every d € N, the dual cone 
Cd(K)* is the convex cone 

(2.6) Arf := |(^^x-#^, aGN3 : <AGM(K)|, 

i.e., the convex cone of vectors ofM.'^^'^^ which have a representing measure 
with support contained in K. 

Proof. For every y = (t/q,) G and / G Crf(K) with coefficient vector 
f G W^'^^: 

(2.7) (y,f) = fc.ya= I ^"^"^"^ = f fdct>> 0. 

Since (12. 7p holds for all / G Crf(K) and all y G A^, then necessarily A^ C 
C7rf(K)* and similarly, C7d(K) C A*. Next, 

A^ = {f GM^W : (f,y) >0 Vy G A,} 

/GM[x]rf : / fd<P>Q V0GM(K) 
^ ASCQ(K), 

and so A^ = Cd(K). Hence the result follows if one proves that A^ is closed, 
because then Crf(K)* = (A^)* = A^, the desired result. So let (y'^) C A^, 
A; G N, with —> y as k ^ oo. Equivalently, Jj^ n.'^dcjjk Ua for all a G N^. 
In particular, the convergence — )• yo implies that the sequence of measures 
{4>k)-, A: G N, is bounded, that is, sup;j(/)fc(K) < M for some M > 0. As K 
is compact, the unit ball of 7V4(K) is sequentially compact in the weak -k 
topology (t(A^(K), C(K)) where C(K) is the space of continuous functions 
on K. Hence there is a finite Borel measure (j) G M(K) and a subsequence 
{ki) such that /j^ gd(t)ki ~^ /k d'^'t' as i ^ oo, for all g G C(K). In particular, 
for every q G NJJ, 

y„ = lim yl = lim yl' = lim / x"#fc^ = / x"#. 
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which shows that y E A^, and so is closed. □ 
And we also have: 

Lemma 2.6. Let K C M" be with nonempty interior. Then the interior of 
Crf(K)* is nonempty. 

Proof. Since Crf(K) is nonempty and closed, by Faraut and Korany [TT| 
Prop. 1.1.4, p. 3] 

int(Q(K)*) = {y : (y, g) > 0, V5 G C,(K) \ {0}}, 

where g € W^'^^ is the coefficient of g G Crf(K), and 

int((:7rf(K)*) / ^ Ca(K) D (-Q(K)) = {0}. 

But g G Crf(K) n (— C(i(K)) implies 5 > and < on K, which in turn 
implies g = because K has nonempty interior. □ 

3. Main result 

Consider the following problem Pq, a non convex generalization of the 
Lowner-John ellipsoid problem: 

Pq.' Let K C be a compact set not necessarily convex and d an even 
integer. Find an homogeneous polynomial g of degree d such that its sublevel 
set Gi := {x : g{'x.) < 1} contains K and has minimum volume among all 
such sublevel sets with this inclusion property. 

In the above problem Pq, the set Gi is symmetric and so when K is a 
symmetric convex body and d = 2, one retrieves the Lowner-John ellipsoid 
problem in the symmetric case. In the next section we will consider the 
more general case where Gi is of the form Gf := {x : g{x — a) < 1} for 
some a G M" and some g G P[x]^. 

Recall that P[x]rf C M[x]rf is the convex cone of nonnegative homogeneous 
polynomials of degree d whose sublevel set Gi = {x : ^(x) < 1} has finite 
volume. Recall also that Crf(K) C IR[x]rf is the convex cone of polynomials 
of degree at most d that are nonnegative on K. 

We next show that solving P is equivalent to solving the convex optimiza- 
tion problem: 



(3.1) V: p= inf 1/ exp(-5) dx : 1 - 5 G Q(K) I . 

Proposition 3.1. Problem P has an optimal solution if and only if problem 
V in iS. 1\) has an optimal solution. Moreover, V is a finite- dimensional 
convex optimization problem. 



Proof. By Lemma [ 

vol(Gi) = , , / exp(-c?)dx 
r(l + n/d) Jj^n 
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whenever Gi has finite Lebesgue volume. Moreover, Gi contains K if and 
only if 1 — (7 G Cd(K), and so P has an optimal solution g* G P[x]rf if and 
only if g* is an optimal solution of V (with value vol(Gi)r(l + n/2{d))). 
Now since g i-> /jj,, exp(— (7)(ix is strictly convex (by Lemma [274|l and both 
C(i(K) and P[x]rf are convex cones, problem P is a finite-dimensional convex 
optimization problem. □ 

We now can state the first result of this paper: Recall that M(K) is the 
convex cone of finite Borel measures on K. 

Theorem 3.2. Let K C M" he compact with nonempty interior and consider 
the convex optimization problem V in \3.1^ . 

(a) V has a unique optimal solution g* € P[x]rf. 

(b) Let g* € P[x]rf be the unique optimal solution ofV. Then there exists 
a finite Borel measure /u* G M(K) such that 

(3.2) f x°exp(-5*)dx = f :>z'^ dfi* , V|a| = d 

(3.3) [{l-g*)dfi* = 0; ^*(K) = ^ / exp{-g*)dx. 

In particular, fi* is supported on the set F := {x G K : g*{x) = 1} (= 
KnGJ) and in fact, ^* can be substituted with another measure v* G Af (K) 
supported on at most C^"*")^""^) contact points ofV. 

(c) Conversely, if g* G ^[x]^; is homogeneous with 1 — g* (z Cd(K.), 
and there exist points (x(i),Ai) G K x M, Aj > 0, i = 1, . . . , s, such that 
g*{x{i)) = 1 for all i = 1, . . . , s, and 



[ x"exp(-g*)dx = V\x(i)°, 



\a\ = d, 



i=l 

then g* is the unique optimal solution of problem V. 



Proof, (a) As ■p is a minimization problem, its feasible set {g G P[x]d : 
1 — g G C(i(K)} can be replaced by the smaller set 

. / exp(-5(x))dx < / exp(-5(o(x))(ix 

g fc -r [XJrf . Jjgn Jjgn 

1 - 5 G Q(K) 

for some g^ G P[x](i. The set F is a closed convex set since the convex 
function g i— )■ Jjg„ exp(— (7)cix is continuous on the interior of its domain. 

Next, let z = (za), a G N^, be a (fixed) element of int(Cd(K)*) (hence 
zo > 0). By Lemma |2 . 61 such an element exists. Then the constraint 1 — 5 G 
Crf(K) implies (z, 1 — 5) > 0, i.e., (z,(7) < zq. On the other hand, being an 
element of P[x]rf, g is nonnegative and in particular g G Crf(K). But then 
by Corollary LI. 6 in Faraut et Koranyi [TT[ p. 4], the set {g G Crf(K) : 
(z, g) < Zq} is compact (where again g G M*^'^) is the coefficient vector of 
g G Cd(K)). Therefore, the set F is a compact convex set. Finally, since 
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g I—)- exp(—g(x))dx is strictly convex, it is continuous on the interior of 
its domain and so it is continuous on F. Hence problem V has a unique 
optimal solution g* G P[x]d. 

(b) We may and will consider any homogeneous polynomial g as an ele- 
ment of M[x](i whose coefficient vector g = {ga) is such that = whenever 
I a I < d. And so Problem V is equivalent to the problem 



(3.4) 



p = 

s.t. 



inf 



exp(-5(x)) dx 

ga = 0, Va e N^; |a| < d 
1-5 € Cd(K), 



where we replaced g S P[x]rf with the equivalent constraints g S M[x]rf and 
5q := for all a € with \a\ < d. Next, doing the change of variable 
h = 1 — g, V reads: 



(3.5) 



V 



P = 
s.t. 



inf 



/ 

Id JR" 



exp(/i(x) — 1) dx 



h. 



a = 0. 
/lo = 1 
h € Crf(K), 



Va G N2; < \a\ < d 



As K is compact, there exists 6 € P[x]rf such that 1 — G int(C(i(K)), i.e., 
Slater's conditioi£| holds for the convex optimization problem V' . Indeed, 
choose X I—)- 0(x) := M~^|[x|['^ for M > sufficiently large so that 1 — 6 > 
on K. Hence with \\g\\i denoting the £i-norm of the coefficient vector of g 
(in M[x]d), there exists e > such that for every h G B{9,e){:= {h G M[x]rf : 
116* — h\\i < e}), the polynomial 1 — /i is (strictly) positive on K. 

Therefore, the unique optimal solution (1— 5*) =: h* G M[x]rf of V in (|3.5p 
satisfies the Karush-Kuhn- Tucker (KKT) optimality condition^, which for 
problem (|3.5p read: 



(3.6) [ x" exp(/i*(x) - l)dx 
Jr" 

(3.7) [ x° exp(/i* (x) -l)dx + 7„ 



(3.8) 



{h*,y*) =0; h*o = 1; K 



y*a, V|al = (i 

Va, y\oi\<d 
0, VO<|al<d 



Slater's condition holds for an optimization problem {min(jo(x) : gk{^) < 0, fc = 
1, . . . , m} if there exists xo such that (7fc(xo) < for all fc = 1, . . . , m. If the gj's are 
convex, continuously differentiable, and Slater's condition holds then the Karush-Kuhn- 
Tucker optimality conditions are necessary and sufficient optimality conditions. 

^For an optimization problem {mmgo{x.) : (/fc(x) < 0, fc = l,...,m}, the KKT- 
optimality conditions hold at a feasible point x* is there exists < A G such that 
V(/(x*) + XifeLi AfcVgfe(x*) = 0, and Xkgki^*) = 0, fc = 1, . . . , m. 
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for some y* = (y*), a G N^J, in the dual cone Cd(K)* C W^'^^ of Crf(K), and 
some vector 7 = (7^), < \a\ < d. By Lemma [231 



Cd(K)* = {y G r ('^) : G M(K) s.t. = [ x" a G }, 

Jk 

and so ()3.2p is just ()3.6p restated in terms of n* . 

Next, the condition {h*,y*) = (or equivalently, (1 — g*,'y*) = 0), reads: 

/ {l-g*)dfi* = 0, 

JK 

which combined with 1 — 5^* G Crf(K) and fi* G M(K), implies that n* is 
supported on K n {x : 5r*(x) = 1} = K n G^. 

Next, let s := Y.\a\=d9*ay*a (= Vo)- From (1 - g*,^*) = 0, the measure 
s~^H* =: V is a probability measure supported on K n G\, and satisfies 
J x.'^dip = s~^y'^ for all |a| = 

Hence by a result of Mulholland and Rogers |22) . there exists an atomic 
measure v* G M(K) supported on K n such that 

/ x"dz^*(x) = / x"dV(x) = s^Sa, V|a| = (i. 

(See e.g. Anastassiou [1, Theorem 2.1.1, p. 39].) As 5*(x) = 1 on the 
support of z/*, we conclude that 



I I , JKnGt JKnGT 



du* 



\a\=d 

i.e., I'* is a probability measure and therefore may be chosen to be supported 
on at most (""'^^"^) points in KnGJ; see [11 Theorem 2.1.1]. Hence in ([32]) 
the measure n* can be substituted with the atomic measure s u* supported 
on at most ("^^~^) contact points in K n G\. 

To obtain ^*(K) = ^ exp(-5-*), muhiply both sides of ([MD-dSZ]) by 
/i* for every a 7^ 0, sum up and use {h* ,y*) = to obtain 



-yl = Y.^lyl = I (r(x)-l)exp(r(x)-l)dx 

5f*(x)exp(-5f*(x))dx 
exp(-5'*(x))dx. 



where we have also used ()2.5p . 

(c) Let /u* := Yld=i where (5x(j) is the Dirac measure at the point 

x(i) G K, i = l,...,s. Next, let y* := / x"d/i* for ah a G N^, so that 
y* G Crf(K)*. In particular y* and g* satisfy 



{i-9*,y*) = I {i-g*)d^i* = 0, 

Jk 
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because g*(-x{i)) = 1 for all i = 1, . . . , s. In other words, the pah' {g*,y*) 
satisfies the KKT-optimality conditions associated with the convex prob- 
lem v. But since Slater's condition holds for V, those conditions are also 
sufficient for g* to be an optimal solution of V, the desired result □ 

Importantly, notice that neither K nor are required to be convex. 

3.1. On the contact points. Theorem 13.21 states that V (hence Pq) has a 
unique optimal solution g* € P[x]rf and one may find contact points x*(«) G 
K n , i = 1, . . . , s, with s < ("^^"^) , such that 

(3.9) = / x"exp(-5*(x))dx = VAax*(i))", \a\ = d, 

for some positive weights X*. In particular, using the identity (|2.5p and 
(1 — g*,y*) = 0, as well as g*{x*{i)) = 1 for all i, 

yo = Y] y*a9*a = 5 / exp(-5r*(x))dx = Vt**- 

\a\=d 1=1 

Next, recall that d is even and let • ~^ M^^"^/^) be the mapping 

X ^ v^/2(x) = (x^*), \a\ = d, 

i.e., the s{d/2)-vector of the canonical basis of monomials (x"^), \a\ = d/2, 
for the homogeneous polynomials of degree d. From ()3.6p . 



(3.10) / V^/2vJ/2exp(-5*)(ix = ^7*V,/2(x*(i))Vrf/2(x*(i))^. 

Hence, when d = 2 and K is symmetric, one retrieves the characterization 
in John's theorem Theorem 2.1], namely that if the euclidean ball ^„ := 
{x : ||x|| < 1} is the unique ellipsoid of minimum volume containing K then 
there are contact points (uj) C ^„ fl K and positive weights (Aj), such that 
AjUjU^ = /„ (where I„ is the n x n identity matrix). Indeed in this 
case, V(^/2(x) = x, (7*(x) = ||x|p and Jj^^ xx"^ exp( — ||x|p)dx = c/„ for some 
constant c. 

So (j3.10p is the analogue for d > 2 of the contact-points property in 
John's theorem and we obtain the following generalization: For d even, let 
l|x||d := (Z]r=i xf)^^^ denote the d-norm with associated unit ball := {x : 
l|x||d < 1}. 

Corollary 3.3. // in Theorem \3.2\ the unique optimal solution G^ is the 
d-unit ball then there are contact points {x*{j)) C K n and positive 
weights Xj, j = 1, . . . , s, with s < ("^^^^); such that for every \a\ = d, 



Y,x*i^*{j)r 



n 

TT / exp(-t'^) if a = 213 



otherwise. 
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4. The general case 

We now consider the more general case where the set Gi is of the form 
{x : g{x - a) < 1} =: Gf where a € M" and g e P[x]rf. 

For every a E M" and g G ^[x]^ (with coefficient vector g E M'*^'^)) define 
the polynomial g^ G M[x]^ by x i-)- (7a(x) := g{x — a) and its sublevel set 
G^ := {x : 5'a(x) < 1}. The polynomial g^ can be written 

(4.1) g^{x) = g{x-a) = ^ pa(a,g)x", 

where g G M*^'^) and the polynomial Pa G M[a, g] is linear in g, for every 
a € N^. 

Consider the following generalization of Pq: 

P; Let K C M" be a compact set not necessarily convex and d G N an 
even integer. Find an homogeneous polynomial g of degree d and a point 
a € M" such that the sublevel set Ga := {x : (^(x — a) < 1} contains K 
and has minimum volume among all such sublevel sets with this inclusion 
property. 

When d = 2 one retrieves the general (non symmetric) Lowner-John el- 
lipsoid problem. For d > 2, an even more general problem would be to find 
a (non homogeneous) polynomial g of degree d such that K C G = {x : 
^(x) < 1} and Gi has minimum volume among all such set G with this 
inclusion property. However when g is not homogeneous we do not have an 
analogue of Lemma 12.21 for the Lebesgue- volume vol(G). 

So in view of ()4.3p . one wishes to solve the optimization problem 

(4.2) V: p= min {vol(G?) : 1 - 5a G ^^(K)}, 

aGR",geP[x]d 

a generalization of (13. ip where a = 0. Let K — a denotes the set {x — a : 
X S K}, and observe that whenever g € P[x]rf, 

(4.3) vol(Gf) = voKG?) = vol(Gi) = ^ [ exp(-5(x)) dx. 

Theorem 4.1. Let K C be compact with nonempty interior and consider 
the optimization problem V in ^.2^ . 

(a) V has an optimal solution (a*,g*) G M" x P[x]rf. 

(b) Let (a*, (7*) G M" x Pfx]^ be an optimal solution ofV. Then there 
exists a finite Borel measure fi* G M(K — a*) such that 



(4.4) / x°exp(-c/*)dx = / x'^ 



dn* , y\a\ = d 



(4.5) 



/ {l-g*)d^i* = 0; ;,*(K-a*) = 5/ exp{-g*)dx. 

JK-a* " JR" 
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In particular, fi* is supported on the set F := {x G K : g*(x — a*) = 1} (= 
KnGf ) and in fact, fi* can be substituted with another measure v* G M(K) 
supported on at most {^^'^~^) contact points ofV. 

Proof. First observe that (j4.2p reads 

(4.6) V: mill I min {vol(Gf) : 1 - 5a 6 Cd{K)]\ , 

and notice that the constraint 1 — i^a £ C'(K) is the same as 1—g G C(K — a). 
And so for every a € M", the inner minimization problem 

min {vol(Gf) : 1 - e Cd{K)} 
sePMd 

of (I4.6P reads 

(4.7) Pa = min {vol(Gi) : 1 - g e Cd{K - a)}. 

sePMd 

From Theorem 13.21 (with K — a in heu of K), problem (14. 7p has a unique 
minimizer g* € P[x]rf with value pa = J^n exp{—g*)dx = J^„ exp( — (7*)dx. 

Therefore, in a minimizing sequence {ai,g^) C x P[x]d, i € N, for 
problem V in (14. 2|) with 

p= lim / exp(-5;)(ix, 

we may and will consider that for every i, the homogeneous polynomial 
5^ S P[x]rf) solves the inner minimization problem (14. TP with a^ fixed. 

As observed in the proof of Theorem 13.21 there is z G int(C(i(K)*) such 
that (1 — g^^jz) > and by Corollary 1.1.6 in Faraut et Koranyi the 
set {h € Cd(K) : (z, /i) < zq} is compact. 

Also, a^ can be chosen with ||a^|| < M for all i (and some M), otherwise 
the constraint 1 — g^i € Cd(K.) would impose a much too large volume 
vol(G^^). 

Therefore, there is a subsequence (ik), A; € N, and a point (a*, 6*) E 
M" X Crf(K) such that 

lim ae^ = a*; lim (5* )^ = 0^, Vq G N^. 

fe— i>oo fe— >oo fc 

Recall the definition (j4.ip of ^^.(x) = (7|(x — a^) for some homogeneous 
polynomial g^ G P[x]rf with coefficient vector g|, i.e., 

(5l*)a = Pa(a|,g|), VaGNS, 

for some polynomials (pa) C M[a, g], a G N^. In particular, for every a G 
with \a\ = d, pa{a},g}) = {gl)a- And so for every a G with \a\ = d, 

K = lim = (5; )„. 
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If we define the homogeneous polynomial g* of degree d by {g*)a = &a 
every a G with \a\ = d, then 

hm (5** )„ = lim po, (a^* , ) , = p« (a* , g* ) , Va G . 

This means that for every a G 

r(x) = (7*(x-a*), xGM". 

In addition, one has the pointwise convergence gl{x) — t- 5'*(x) for all x G M". 
Therefore, by Fatou's Lemma (see e.g. Ash jt2j), 

p = lim / exp(— g'^)dx > / liminf exp(— (/|) dx = / exp(— dx, 

JK" JR" JR" 

which proves that (a* ,g*) is an optimal solution of ()4.2p . 

In addition g'* G P[x]rf is an optimal solution of the inner minimization 
problem in (I4.7p with a := a*. Otherwise an optimal solution h G P[x](i 
of (j4.7p with a = a* would yield a solution (a* , h) with associated cost 
J^n exp(— /i) strictly smaller than p, a contradiction. 

Hence by Theorem 13.21 (applied to problem (14. 7p ) there is a finite Borel 
measure p* G M(K — a*) such that 

/ x" exp(— 5'*)(ix = / x." dp, V|a| = d 

7ri Jk-sl* 

[ {l-g*)dp* = 0; MK-a*) = 5/ exp(-5*)a!x. 

And so ;U* is supported on the set 

Fa* = {x G K - a* : 5*(x) = 1} = {x G K : r7*(x - a*) = 1}. 

Invoking again [Tl Theorem 2.1.1, p. 39], there exists an atomic measure 
1^* G M(K - a*) supported on K - a* n such □ 



5. A COMPUTATIONAL PROCEDURE 

Even though V in (j3.ip is a finite-dimensional convex optimization prob- 
lem, it is hard to solve for mainly two reasons: 

• From Theorem 12.41 the gradient and Hessian of the (strictly) convex 
objective function g J exp{—g) requires evaluating integrals of 
the form 

/ x" exp(-5(x)) dx, Va G N^, 

JR" 

a difficult and challenging problem. (And with a = one obtains 
the value of the objective function.) 

• The convex cone (7^(K) has no exact and tractable representation to 
efficiently handle the constraint 1 — g £ Crf(K) in an algorithm for 
solving problem (|3.ip . 
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However, below we outline a numerical scheme to approximate to any desired 
e-accuracy (with e > 0): 

- the optimal value p of (13. ip . 

- the unique optimal solution g* G P[x]d of V obtained in Theorem 13.21 

5.1. Concerning gradient and Hessian evaluation. To approximate 
the gradient and Hessian of the objective function we will use the following 
result: 

Lemma 5.1. Let g G P[x]rf and let Gi = {x : ^(x) < 1}. Then for every 
a G N" 

(5.1) / x"exp(-g)rfx = r fl+ ""^J"* ^ / x" dx. 

The proof being identical to that of Lemma 12.21 is omitted. So Lemma 
15.11 relates in a very simple and explicit manner all moments of the Borel 
measure with density exp(— 5^) on R" with those of the Lebesgue measure 
on the sublevel set Gi. 

It turns out that in Henrion et al. |14] we have provided a hierarchy of 
semidefinite program^ a convex optimization problem to approximate as 
closely as desired, any finite moment sequence {za)-, a G N", defined by 

Zc, = / x" dx, Q G N^. 

where O is a compact basic semi-algebraic set of the form {x : 0j(x) > 0, j = 
1, . . . , m\ for some polynomials {9j) C M[x]. For more details the interested 
reader is referred to [14j. 

Hence in any minimization algorithm for solving V, and given a current 
iterate g G P[x]rf, one may approximate as closely as desired the value at 
g of the objective function as well as its gradient and Hessian by using the 
methodology described in [HI [19]. 

5.2. Concerning the convex cone (7rf(K). We here assume that the com- 
pact (and non necessarily convex) set K C is a basic semi-algebraic set 
defined by 

(5.2) K = {x G M!^ : u;,(x) > 0, j = 1, . . . , s}, 

for some given polynomials {wj) C M[x]. Denote by C M[x]2fc the convex 
cone of SOS (sum of squares) polynomials of degree at most 2k, and let wq 
be the constant polynomial equal to 1, and Vj := [deg(u;j)/2] , j = 0, . . . , s. 



A semidefinite program is a finite-dimensional convex optimization problem which in 
canonical form reads: minx{c-^x : Aq + X^^^i A.kXk ^ 0}, where c £ R', the A^'s are 
real symmetric matrices, and the notation A ^ stands for A is positive semidefinite. 
Importantly, up to arbitrary fixed precision it can be solved in time polynomial in the 
input size of the problem. 
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With k fixed, arbitrary, we now replace the condition 1 — (7 € Cd(K.) with 
the stronger condition 1 — g & Ck {c Crf(K)) where 

f s 

(5.3) '^'^ = ) '^j ■ ^ ^k~vj , i = 0, 1, . . . , s 

It turns out that membership in Ck translates into Linear Matrix Inequali- 
ties (LMIs) on the coefficients of the polynomials g and cr^'s; see e.g. [19J. 
If K has nonempty interior then the convex cone Ck is closed. 

Assumption 1 (Archimedean assumption). There exist M > and /c E N 
such that the quadratic polynomial x 1— ?> Af — ||x|p belongs to C^- 

Notice that Assumption [T] is not restrictive. Indeed, K being compact, 
if one knows an explicit value M > such that K C {x : ||x|| < M}, then 
its suffices to add to the definition of K the redundant quadratic constraint 

Ws+i{^) > 0, where Ws+i(x) := — ||x|p. 

00 

Under Assumption [H Crf(K) = \^ Ck, that is, the family of convex cones 

fc=0 

(Ck), A; G N, provide a converging sequence of (nested) inner approximations 
of the larger convex cone C(i(K.)- 

5.3. A numerical scheme. In view of the above it is natural to consider 
the following hierarchy of convex optimization problems (Vk), A: € N, where 
for each fixed k: 



Pk = min / exp(-5f) dx 



(5.4) s.t. l-g = ^ajWj 

j=0 

ga = 0, V|a| < d 

9 e I^Md; aj G ^k-vp J = 0, . . . , s. 

Of course the sequence (pk), A; G N, is monotone non increasing and pk ^ P 
for all k. Moreover, for each fixed A G N, 7^^ is a convex optimization 
problem which consists of minimizing a strictly convex function under LMI 
constraints. 

From Corollary 12.31 J^n exp{—g)dx < 00 if and only if 5 G P[x]d and 
so the objective function also acts as a barrier for the convex cone P[x]rf. 
Therefore, to solve Vk one may use first-order or second-order (local min- 
imization) algorithms, starting from an initial guess go G P[x]rf. At any 



Linear Matrix Inequality (LMI) is a constraint of the form Ao + X]^=i ^i^e h 
where each A^, £ = 0, . . . ,t, is a, real symmetric matrix and the notation A ^ stands for 
A is positive semidefinite. Equivalently A(x) >^ where A(x) ~ Ao + X]£=i ^e^e, a real 
symmetric matrix whose each entry is afline in x G R*. An LMI always define a convex 
set, i.e., the set {x G R* : A(x) ^ Ao} is convex. 
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current iterate g G P[x]d of such an algorithm one may use the method- 
ology described in [14] to approximate the objective function J exp(— (7) as 
well as its gradient and Hessian. 

Theorem 5.2. Let K in ^5.2\) be compact with nonempty interior and let 
Assumption[I\hold. Then there exists ko such that for every k > k^, problem 
Vk in |5.^[ j has a unique optimal solution G P[x](i. 

Proof. Firstly, Vk has a feasible solution for sufficiently large k. Indeed 
consider the polynomial x go{'^) = Y17=i^i which belongs to P[x]rf. 
Then as K is compact, M — go > on K for some M and so by Putinar's 
Positivstellensatz [26J, 1 — go/M G Ck for some ko (and hence for all /c > /cq)- 
Hence go/M is a feasible solution for Vk for all k > k^. Of course, as 
Ck C Cd(K), every feasible solution g G P[x]rf satisfies < g < 1 on K. 
So proceeding as in the proof of Theorem 13.21 and using the fact that Ck is 
closed, the set 

{5 G P[x]d nCfc : / exp{-g)dx< exp(-|^)(ix}, 

Jr" Jm." 

is compact. And as the objective function is strictly convex, the optimal 
solution (7^ G P[x]rf n Ck is unique (but the representation of 1 — g^ in (15. 4p 
is not unique in general). □ 

Remark 5.3. // desired one may also impose g to be convex by simply 
requiring z^V^5'(x)z > for all (x, z). 

We now consider the asymptotic behavior of the solution of (15. 4p as A; — ?> 
00. 

Theorem 5.4. Let K in i5.2\) be compact with nonempty interior and let 
Assumption m hold. If p (resp. pk) is the optimal value of V (resp. Vk) 
then p = lim pk- Moreover, for every k > k^, let g^ G P[x]rf be the unique 

k^oo 

optimal solution of Vk- Then as k ^ 00, g^ — ?• g* where g* is the unique 
optimal solution ofV. 



Proof. By Theorem 13.21 V has a unique optimal solution g* G P[x](;;. Let 
e > be fixed, arbitrary. As 1 — g* € Cd(K.), the polynomial 1 — g/{l + e) is 
strictly positive on K, and so by Putinar's Positivstellensatz [26], l — g*/{l + 
e) belongs to Ck for all k > k^ for some integer k^. Hence the polynomial 
5*/(l + e) G P[x]d is a feasible solution of Vk for all k > k^. Moreover, by 
homogeneity, 

/ exp(--^)(ix = {l + eTl" f exp(-5*)dx 
= (l + e^/V 
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This shows that /Ofc < (1 + e)"'^'^p for ah k > k^. Combining this with pk > p 
and the fact that e > was arbitrary, yields the convergence pk ^ p as 
k — )• oo. 

Next, let y G int(Crf(K)*) be as in the proof of Theorem 13.21 From 
1 — £ Ck we also obtain (y, 1 — g^) > 0, i.e.. 

Recall that the set {g G Cd(K) : (y,^) < yo} is compact. Therefore there 
exists a subsequence (ki), ^ € N, and g E Crf(K) such that g^^ g as 
£ — > oo. In particular, 1 — g £ Crf(K) and ga = whenever |a| < d (i.e., g is 
homogeneous of degree d). Moreover, one also has the pointwise convergence 
lim£_!.oo (7^^ (x) = ^(x) for all x G M". Hence by Fatou's lemma, 

p = lim Pk, = lim / exp(-5^ (x))(ix 
> / liminfexp(-fi(^ (x))(ix 

= / exp(-^(x))(ix = p, 

which proves that g is an optimal solution of V, and by uniqueness of the 
optimal solution, g = g*. As {gkf), ^ € N, was an arbitrary converging 
subsequence, the whole sequence {gD converges to g*. □ 

6. Conclusion 

We have considered non convex generalizations Pq and P of the Lowner- 
John ellipsoid problem where we now look for an homogeneous polynomial 
g of (even) degree d > 2. Importantly, neither K not the sublevel set G 
associated with g are required to be convex. However both Pq and P have an 
optimal solution (unique for Pq) and a characterization in terms of contact 
points in K n G is also obtained as in the Lowner-John ellipsoid theorem. 
Crucial is the fact that the Lebesgue volume G is a strictly convex function 
of the coefficients of g. This latter fact also permits to define a hierarchy 
of convex optimization problems to approximate as closely as desired the 
optimal solution of Pq. 
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